Resonance Region Structure Functions and Parity Violating Deep Inelastic Scattering 
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The primary motive of parity violating deep inelastic scattering experiments has been to test the 
standard model, particularly the axial couplings to the quarks, in the scaling region. The measure- 
ments can also test for the validity of models for the off-diagonal structure functions F^ z (x,Q 2 ) 
in the resonance region. The off-diagonal structure functions are important for the accurate cal- 
culation of the 7^-box correction to the weak charge of the proton. Currently, with no data to 
determine F^f 3 (x, Q 2 ) directly, models are constructed by modifying existing fits to electromag- 
netic data. We present the asymmetry value for deuteron and proton target predicted by several 
different _F7f 3 (x, Q 2 ) models, and demonstrate that there are notable disagreements. 



I. INTRODUCTION 

Parity violating deep inelastic scattering experiments 
(PVDIS) in the past Q and still today have been mo- 
tivated by the desire to search for physics beyond the 
Standard Model in lepton-quark neutral current inter- 
actions. Consistent with the focus on interpreting the 
results in terms of quark couplings, the kinematics are 
chosen to mainly lie in the scaling region. 

However, the recent interest in the larger than ex- 
pected 7-Z box correction to the Qweak experiment Q 
provides a further motivation for PVDIS, particularly for 
any data that may lie in the resonance region. QWeak is 
also an experiment designed to test the standard model, 
using elastic electron-proton scattering with polarized 
electrons, and obtaining an accurate result requires good 
knowledge of higher order corrections. Hence the interest 
in the 7-Z boxes. 

Both PVDIS and the 7-Z box calculations are depen- 
dent upon the off-diagonal structure functions ^2 3> de- 
fined from the spin-averaged tensor 
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Functions F^ Z come from the vector terms in the Z- 

boson current and F^ z comes from the axial vector 
terms; we may call them Zy and Za couplings for brevity. 
The 7-Z box, Fig. [TJ can be calculated dispersively in 

The surpris- 



(1) 



terms of the structure functions F^ z 3 
ing recent result Q was that the vector Z-boson contri- 
butions were larger than expected and numerically com- 
parable to the axial contributions, which had comprised 
the bulk of the earlier estimates (see for example 0). 
In the scaling region and with the standard model one 
obtains the 7Z structure functions from parton distribu- 
tion functions. However, the contributions from F± z are 
given as integrals that have their main support in the 



resonance region and at moderate Q 2 . In the absence of 
data for these structure functions, one has obtained them 
using model-based modification of data from other chan- 
nels, in particular starting from the standard electromag- 
netic structure functions, here called F^ 2 - Regarding the 

F^ z integrals, the main support is at higher energy and 
higher Q 2 , so that known parton distribution functions 
can be used to obtain the bulk of the axial contributions. 
However, one still wants to know F^ z in the resonance 
region. In principle, this can be obtained from the charge 
current reaction data, since only the Weak axial current 
is involved. However, Weak interaction resonance region 
data is scarce so that modeling is still needed. 

PVDIS allows a direct measurement of the 7Z struc- 
ture functions. The PVDIS asymmetry is given by fic| 
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Here, A is the incoming lepton helicity, E and E' are 
the incoming and outgoing lepton energy in the target 





FIG. 1: The 7-Z box diagrams. 
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rest frame, v = E — E', 9 is the lcpton scattering angle, 
x = Q 2 / (2Mv), y = v/E, and in the standard model, 
g e A = —1/2 and g v = — 1/2 + 2 sin 2 9w, where 6w is the 
Weinberg angle. 

In principle, measuring the asymmetry over a range of 
angles, initial energies, and Q 2 allows a full determina- 
tion of the structure functions. In practice, at least at the 
outset, the data will be limited, and so it will be useful 
to predict the asymmetries in the resonance region using 
the models developed for evaluating the 7Z boxes, and 



then use the data to spot-check the models. It is worth 
mentioning at the outset that for most relevant kinemat- 
ics the numerical contribution of the F^ z term will be 
small. 



Although this article is focused on the resonance re- 
gion, for the purpose of contrast and comment we give 
the PVDIS asymmetry formula specialized to the scaling 
region, 



Apvdis - 3GfQ2 2Clu ( UA + ~ C ^(d A + d A + s A + s A ) + Y (2C 2u (ua - u A ) - C 2d (d A - d A )) 
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where q A is the distribution function for quark q in target 
A. The F^ z term has become the term containing Y(y) 
and C 2q , where 
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In the standard model, gy 
W, and g\ = § 
The F^ z contribution is largest for x — > 1 and y 1, 
where one has Y 1, and expects the antiquark and 
strange quark distributions to be very small. For the 
standard model in this limit, one expects the C 2q terms 
to be about 12% of the C\ q terms in the deuteron, where 
one can let u A = d A . For the proton, the effect is some- 
what larger but dependent on the down to up quark ra- 
tio in the valence region (x —¥ 1). Beyond the standard 
model searches will work in the scaling region and seek 
deviations from this result. 

We, on the other hand, will here accept the standard 
model, and work in the resonance region and hope to 
learn about F^ z . 

We continue by showing predictions for the PVDIS 
asymmetry in the resonance region based on several mod- 
els that have been proposed and used in the box cal- 
culations, and discussing the reasons for the differences 
among these model predictions, and then offer some con- 
clusions. 



II. ASYMMETRIES IN THE RESONANCE 
REGION 

Figures [2] and [3] show several examples of what models 
predict for the PVDIS asymmetry in the resonance re- 
gion, choosing for definiteness the two Q 2 values where 
the 6 GeV PVDIS experiment has been run. 



The models for the jZ structure functions in the res- 
onance region have been mainly discussed in the context 
of a proton target, and so we begin with the proton in 
Fig. [2J The top two panels show results for the PVDIS 
asymmetry that follow from four different models of how 
to convert the electromagnetic structure functions to the 
j-Z ones. The vertical dashed line in each figure shows 
the value of W which is targeted in the current (deuteron) 
experiment. 

Each of the jZ structure function models that we show 
have their roots in the Christy-Bosted fit [ll| to 77 struc- 
ture function data. To convert to the j-Z structure func- 
tions, each of the fit's seven electromagnetic resonance 
contributions is modified by a corrective ratio. Explicitly, 
for the transverse contributions, one forms the ratios 



Cr 
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(where A\ is the transverse helicity amplitude and A is 
the helicity magnitude of the resonance), and multiplies 
each respective electroproduction resonance contribution 
by its Cr(Q 2 ). The jZ structure function models differ 
in the way they obtain the corrective ratios, as well as in 
their handling of the nonrcsonant background. 

In the top panels, the curve labeled "CQM" follows 
from the model used by us in 0. The model is a con- 
stituent quark model for resonance production and is 
used to calculate helicity amplitudes for generic vector 
couplings. Helicity amplitudes for Z-boson (photon) ex- 
change are obtained by inserting the g v vector Weak 
charges [e q quark charges). Many transitions receive con- 
tributions from both electric and magnetic excitations. 
With extra information about the helicity behavior of 
the spin- 3/ 2 resonances, one can determine the individ- 
ual contributions from the electromagnetic data and then 
convert to the Zy amplitudes. The background we use 
comes from the background part of the Christy-Bosted 
fits to inelastic electron scattering data, modified for the 
^Z case. 
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FIG. 2: Normalized proton asymmetry for Q 2 — 1.1 GeV 2 and 1.9 GeV 2 as a function of W. The top panels show results for 
several parameterizations, one based on a constituent quark model based modification of the Christy-Bosted electromagnetic fits 
(black solid), one obtained using MAID fits to the resonance helicity amplitudes (green dotted), and two suggestions following 
Gorchtein et al. (red dashed and blue dash-dot). The middle and bottom panels give uncertainty limits for the constituent 
quark model and MAID cases, respectively, with the grey band being uncertainty for the nonresonant terms alone, and the pink 
band adding uncertainty in the resonant contributions to obtain the total. The dashed vertical lines indicate the kinematic 
points for the 6 GeV PVDIS (deuteron) experiment; each corresponds toifs 0.3. 



The curves labeled "Model I" and "Model II" both 
follow from the models used in (GHRM). The analy- 
sis of the corrective ratio, Eq. ((6j), for GHRM begins by 
considering the relation 

2 (R+ 1 Jl - \p) = ( 1 - 4 sin 2 e w ){R+\j; \p) 

- (R°\J^\n) ~ ( R+ MP) ■ ( ? ) 

The strange quark contribution is generally neglected, as 
suggested by data fl2T - fl4l ]. so that the vector Z-current 
helicity amplitudes can be gotten from neutron and pro- 
ton data. GHRM implement this only at the photopro- 
duction point, using proton and neutron helicity ampli- 
tudes given by the Particle Data Group They as- 



sume that the Q 2 dependence of the Zy and electromag- 
netic matrix elements are the same, so that they have Q 2 
independent Cr's. They obtain resonant electromagnetic 
matrix elements from the resonant part of the Christy- 
Bosted fit [llj|, slightly modified to better match their 
background choices. 

The GHRM backgrounds come from two fits to electro- 
production in a higher energy diffractive region, pushed 
to lower energy and isospin modified for the jZ case. 
Model I is based on a color dipole model based fit of 
Cvetic et al. [l5j]. In this model, the photon fluctuates 
into a q-q pair which interacts with the proton via gluon 
exchange. Cvetic et al. obtain a functional form with pa- 
rameters, which are constrained by data, and a good fit 
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FIG. 3: Deuteron asymmetry for Q 2 = 1.1 GeV 2 and 1.9 GeV 2 as a function of W c a, where W 2 H = M 2 + 2Mv - Q 2 . The top 
panels show results for several parameterizations, one based on a constituent quark model based modification of the Christy- 
Bosted electromagnetic fits (black solid), one obtained using MAID fits to the resonance helicity amplitudes (green dotted), and 
one using photoproduction data in the style of Gorchtein et al. (red dashed). The middle and bottom panels give uncertainty 
limits for the constituent quark model and MAID cases, respectively, with the grey band being uncertainty for the nonresonant 
terms alone, and the pink band adding uncertainty in the resonant contributions to obtain the total. The dashed vertical lines 
again indicate the kinematic points for the 6 GeV PVDIS experiment; each corresponds to x ~ 0.3. 



is obtained for low Q 2 , high energy structure functions. 
GHRM 0) extrapolate this to model the background in 
the resonance region, and obtain the j-Z structure func- 
tions by changing the averaged photon-quark coupling to 
Z-quark coupling. Model II is based on a generalized vec- 
tor meson model based fit of Alwall and Ingelman fl6| . 
This model fits the electromagnetic structure functions 
at low Q 2 and high energy by relating them to the to- 
tal photon proton cross sections and coupling the photon 
through vector meson intermediaries, given in terms of p, 
w, and (f> plus a background or continuum contribution. 
Similar to their Model I modifications, GHRM extrap- 
olate the Model II fit to lower energies and use it as a 
background in the resonance region, first showing ade- 



quate fits to purely electromagnetic data and then mod- 
ifying the fit for the 7-Z case by considering Z-vector 
meson as well as 7- vector meson mixings. 

The curve labeled "MAID" used the MAID fits to 
the resonance electroproduction helicity amplitudes with 
proton and neutron targets fl7l ] to obtain the vector Z- 
boson matrix elements from Eq. ([7]) and thence a fully 
Q 2 dependent corrective ratio. These are then used to 
transform each of the seven resonant contributions in the 
Christy-Bosted fit into the ^Z structure function terms. 
The background is obtained as previously described when 
using the CQM model. 

The F$ z terms have small numerical impact here, at 
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most 3% of the total. The F^ z contributions are not big 
even in the scaling region under optimal kinematics, and 
at kinematics relevant to this paper, the function Y(y) 
is well below unity and the presence of sea quark and 
strange quark contributions further reduce the relative 
size of the F^ z terms. To justify the numerical result, 
we need some estimate for F^ z . For the resonance contri- 
bution we use a corrective ratio with a helicity amplitude 
for the Z-boson axial coupling and once again evaluate 
the amplitudes with our constituent quark model. We 
have also obtained resonance contributions from the four- 
resonance fits of [HI], and find the difference it makes in 
our plots is very slight. For the background one takes 
guidance from the scaling region. For very low-x, all 
light quarks and antiquarks have similar distributions, 
so that F^ z , which depends on differences q(x) ~ q(x), 
is about zero compared to F^ 1 . Alternatively, in a va- 
lence quark dominated region where the up quark dis- 
tribution is twice the down quark distribution, one has 
F~f IF^ 1 = 10/3. We take the average of the two limits, 
taking F^ z for the proton to be 5/3 the Christy-Bosted 
-F 77 (with 100% uncertainty bounds on this term, when 
we discuss uncertainty bounds below, to accommodate 
the two limits). 

One can see the results of the various models differ. 
This has much to do with the Q 2 dependences of the 
resonance amplitudes in the different models. Note the 
deep dip in the quark model in the second resonance re- 
gion contrasting with the deep dip in the third resonance 
region for the GHRM predicted asymmetry. GHRM sets 
the ratio of the vector Z to electromagnetic matrix ele- 
ment by their values at Q 2 — 0. The quark model consid- 
ers that both the second and third resonance regions have 
a spin-3/2 or higher resonance, which each have two he- 
licity amplitudes and these amplitudes have different Q 2 
falloffs. The weightings of the two amplitudes changes 
when switching from the photon to the Z-boson case, 
which leads to a different overall falloff with Q 2 for the 
Z-boson contributions, faster for the second resonance 
region and slower for the third resonance region. 

The deep dip seen in the MAID based asymmetry is 
due to the behavior of the Roper resonance in the MAID 
model. In the MAID fits, the proton transverse helic- 
ity amplitude changes sign, at about 2/ 3 GeV 2 , while the 
neutron does not, leading to a sign change in some of 
the interference terms giving the 7Z structure functions 
following from Eq. The negative contribution from 
the Roper leads to the dip. The GHRM and the quark 
model (as it happens for the Roper) have a Q 2 indepen- 
dent electromagnetic to Z boson amplitude and hence no 
sign change and no dip in the Roper region. 

For the uncertainty limits of the smooth non-resonant 
background to i*\ 7 f m the quark model, we take guid- 
ance from scattering off collections of quarks with scant 
final state interactions. In a full SUf(3) limit, where all 
light quarks are equally likely and which may be perti- 
nent in a high-energy x = Q 2 /(2Mv) —¥ limit, one has 



F±2 /Fi2 = 1 + Qw ■ hi a valence quark limit with 
SU(6) wave functions, one gets (2/3 + Q^ LO ) for the 
same ratio. The latter is better at high- a; and the former 
is better at low- a; and we take the mean for our central 
curve, and use the extremes to set the uncertainty esti- 
mates which are shown in the four lower panels of Fig [5] 
Regarding the conversion of the resonances from electro- 
magnetic to Z-boson matrix elements, we assigned a 10% 
uncertainty in each matrix element [6|. The uncertain- 
ties are shown in the middle panels of Fig. [2] are for the 
CQM, with the grey band showing the uncertainty limits 
for the nonresonant terms, and the pink band adding un- 
certainty from the resonant contributions to obtain the 
total. The bottom panels show similar bands for the 
MAID based asymmetry predictions; results for GHRM 
are qualitatively similar. 

For the deuteron asymmetry, Fig. [31 wc modify the 
Bosted-Christy deuteron fits to electromagnetic structure 
function data (l9j . Extending the quark model to obtain 
the Z current matrix elements on a neutron target is of 
course straightforward. Also, the Weak isospin rotation 
giving the Z boson on neutron matrix elements reflects 
the proton case, 

2(R°\J z -\n) Asm 2 8 w )(R°\j;\n) 

-(R + \j;\p)-(R°\si,s\n). (8) 

The constituent quark model and MAID based treat- 
ments of the resonances are thus both straightforward to 
extend to the deuteron. However, a problem in extend- 
ing the GHRM treatment to the deuteron is that we have 
no equivalent of the color dipole model [l5| or general- 
ized vector dominance fl6j fits, which were used for the 
background, to the deuteron. However, we can, as we do 
for the quark model and MAID, use the Bosted-Christy 
deuteron background, suitably modified for the 7Z case, 
and combine with the GHRM resonance fits to do a fit 
in their style. 

The F^ z structure function is again a small contribu- 
tion, and we obtain an estimate of it in the same way 
as we did for the proton target. The background part 
of F^ Z I F^ 1 is still essentially zero in high-energy low-x 
limit, but becomes 18/5 for the deuteron in the valence 
only limit. We again take the average of the two limits, 
taking F^ z for the proton to be 9/5 the Bosted-Christy 
F 77 for the deuteron, with 100% uncertainty bounds to 
accommodate the two limits. 

Comments about the reasons for the differences be- 
tween the quark, GHRM, and MAID based treatments 
are similar to those made for the proton case. 

III. CONCLUSION 

Parity violating deep inelastic asymmetry measure- 
ments were proposed as a source of interesting and useful 
information about the 7Z interference structure func- 
tions in the scaling region. The intent was to use a 
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deuteron target to minimize the uncertainty in remov- 
ing effects of the large F^ z and F^ Z terms and hence 
isolate F% . In the latter, the -Z-quark coupling is axial 
vector, and the goal was to measure if the coupling was 
in accord with the standard model. 

On the other hand, if the standard model is valid, 
PVDIS provides a way to measure the structure func- 
tions in any kinematic region. In particular, the results 
in the resonance region would be very interesting and 
useful. The interest in this region is to measure or con- 
strain the F^ z structure functions, and the smallness of 
the F^ z becomes an advantage. The resonance region 
F±2 are n °t in fact well predicted from existing data in 
other channels. Several models are available, and we have 
shown the PVDIS asymmetries that follow from several 
of them, and note that the predicted asymmetries are 
fairly distinct, differing by several tens of percent in some 
kinematic regions. 

The 7.Z structure functions have an immediate appli- 
cation, which is in the calculation of the box correc- 
tions to elastic parity violating electron-proton scatter- 



ing [4|-l7| . The corrections are calculated using dispersion 
relations, and the numerical contributions from the Z 
boson having vector or axial interactions is about the 
same. However, the weightings within the integrals are 
such that when the Z has vector couplings, which give the 
F^ z , the main support comes from the resonance region. 
Hence the direct interest in measuring these functions 
there. For the Z-boson axial couplings, which give F^ z , 
most of the "/Z box contribution comes from the scaling 
region Q . Hence there is less need from the present view- 
point for measuring F^ z in the resonance region, and the 
smallness of the F^ z contribution to PVDIS can even be 
counted as an advantage. 
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